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ABSTRACT 



> ' 

in : 

|~~» Aims. The development of a general code for 3D Lya radiation transfer in galaxies to understand the diversity of Lya line profiles observed in 
f^") ■ star forming galaxies and related objects. 

OO ' Methods. Using a Monte Carlo technique we have developed a 3D Lya- radiation transfer code allowing for prescribed arbitrary hydrogen 
density, ionisation, temperature structures, and dust distributions, and arbitrary velocity fields and UV photon sources. 

Results. As a first test and application we have examined the Lya line profiles predicted for several simple geometrical configurations and 
their dependence on the main input parameters. Overall, we find line profiles reaching from doubly peaked symmetric emission to symmetric 
Voigt (absorption) in static configurations with increasing dust content, and asymmetric red- (blue-) shifted emission lines with a blue (red) 
O |' counterpart ranging from absorption to emission (with increasing line/continuum strength) in expanding (infalling) media. In particular we 
find the following results which are of interest for the interpretation of Lya profiles from galaxies. 1 ) Standard Lya absorption line fitting of 
global spectra of galaxies may lead to an underestimate of the true hydrogen column density in certain geometrical conditions. 2) Normal 
(inverted) P-Cygni like Lya profiles can be obtained in expanding (infalling) media from objects without any intrinsic Lya emission, as a 
^ natural consequence of radiation transfer effects. 3) The formation and the detailed shape of Lya profiles resulting from expanding shells has 
■ been thoroughly revised. In particular we find that, for sufficiently large column densities (7V H ^ 10 20 cm -2 ), the position of the main Lya 
. m emission peak is quite generally redshifted by approximately twice the expansion velocity. This is in excellent agreement with the observations 
, of z ~ 3 LBGs, which show that Lya is redshifted by ~ 2V* cxp , where V* cxp is the expansion velocity measured from the interstellar absorption 
» I ■ lines blueshifted with respect to the stellar redshift. This finding indicates also that large scale, fairly symmetric shell structures must be a good 
' description for the outflows in LBGs. 
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1. Introduction 

The Lya line plays an important role in a variety of astro- 
physical problems, especially as a diagnostic tool to observe 
and study the high redshift universe. It is a simple redshift in- 
dicator for distant galaxies, a frequently used star formation 
rate diagnostic at high z, as well as an important tool probing 
the ionisation state of the intergalactic medium and hence the 
reionisation epoch. The Lya emission line is e.g. a strong fea- 
ture observed in nearby star forming galaxies, distant Lyman 
break galaxies (LBGs), sub-mm galaxies, emission line se- 
lected galaxies (LAE, for Lyman-a emitters), and in the enig- 
matic so called Lyman-a blobs (LABs) whose nature remains 
debated (e.g. Steidel et al. 2000, Dijkstra et al. 2005b). 

Since the early suggestion of strong Lya emission from 
young high redshift galaxies by Partridge & Peebles (1967) and 
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until the late 1990s, only few Lya emitters have been found (cf. 
Djorgovski & Thompson 1992). This lack of Lya emission de- 
tection has triggered a variety of studies discussing the possi- 
ble physical effects (mostly metallicity, dust, neutral gas kine- 
matics, and geometry) which may significantly affect and sup- 
press the Lya emission and the resonance line radiation trans- 
fer, thereby reducing the observed Lya intensity and destroying 
simple expected correlations, e.g. between Lya intensity and 
metallicity, Lya intensity and UV continuum flux and others 
(Meier & Terlevich 1981, Hartmann et al. 1988, Neufeld 1990, 
Chariot & Fall 1993, Valls-Gabaud 1993, Kunth et al. 1998, 
Tenorio-Tagle et al. 1999, Mas-Hesse et al. 2003). 

In the last few years, with the availability of deeper and 
wider surveys such as the Large Area Lyman Alpha (LALA) 
survey and the Subaru Deep Field survey, many emission 
galaxies have been detected (cf. Hu et al. 1998, 2004; Kudritzki 
et al. 2000, Rhoads et al. 2000, Ouchi et al. 2003, Taniguchi 
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et al. 2005). Although the majority of these distant Lya emit- 
ters shows rather simple asymmetric line profiles, the overall 
diversity of the observed Lya line shapes, both from star form- 
ing galaxies in the nearby universe and at high-z, is quite het- 
erogeneous and complex. The observed line profiles include 
schematically pure Voigt absorption profiles, P-Cygni profiles, 
double peak profiles, pure (symmetric) emission line profiles, 
and combinations thereof (see e.g. Kunth et al. 1998, Mas- 
Hesse et al. 2003, Shapley et al. 2003, Moller et al. 2004, 
Venemans et al. 2005, Wilman et al. 2005, Noll et al. 2004, 
Tapken et al. 2004, Tapken 2005). 

Although in principle the main physical processes shaping 
the Lya line are known, in practice the inferences drawn so 
far from Lya observations rely mostly on rather simple mea- 
surements (e.g. line flux) or on oversimplified Voigt-profile fits, 
which often have no strong physical motivations. For example, 
for a Lya line profile formed purely within a galaxy (i.e. ne- 
glecting subsequent alterations from the intergalactic medium 
and/or intervening clouds), it is physically inconsistent to fit 
one or several Voigt profiles without making strong implicit as- 
sumptions on the geometry of the neutral gas. In such a case, 
due to the resonance scattering nature of Lya it is in fact un- 
likely that the resulting emergent line profile is actually a Voigt 
profile and a detailed radiation transfer calculation has thus to 
be carried out to predict the proper shape of the emergent res- 
onance line profile. In general, quantitative simulations for ap- 
propriate geometries and gas kinematics, which take properly 
into account the main physical processes of Lya line formation 
and radiation transfer are therefore needed for a better under- 
standing of the variety of observed Lya line profile morpholo- 
gies. The physical properties of the Lya emission mechanism 
and the ones of the ambient gas hosting the Lya emitter, which 
both shape the observed Lya line, must be properly modeled to 
investigate the correspondent impact on the emergent line pro- 
file and their possible degeneracies. This in turn will provide 
useful hints to guide the interpretation of the observed profiles 
and to gain insight on the physical properties of the Lya emit- 
ters and of their environment. 

Analytic solutions for the Lya radiation transfer problem 
have been derived for simple geometries. Neufeld (1990) has 
extensively studied the case of static plane parallel slabs, yield- 
ing important insight on the line formation mechanism and pro- 
viding solutions for configuration including dust, Bowen fluo- 
rescence and Lya pumping of H2 Lyman band lines. The case 
of a static, uniform sphere has been recently studied by Dijkstra 
et al. (2005a). Loeb & Rybicky (1999) and Rybicky & Loeb 
(1999) have derived solutions for Lya scattering in a Hubble 
flow. However, more general geometries and velocity fields do 
not allow for an analytic solution and require alternative ap- 
proaches. 

Over the last few years, several groups have developed new 
numerical algorithms, mostly based on Monte Carlo techniques 
(Spaans 1996, Ahn et al. 2001, 2002; 2003, Ahn 2004; Zheng 
& Miralda-Escude 2002; Richling et al. 2001; Richling 2003; 
Kobayashi & Kamaya 2004; Cantalupo et al. 2005; Dijkstra 
et al. 2005ab; Hansen & Oh 2006; Tasitsiomi 2006). Some of 
these codes have been specifically designed for and can be re- 
liably applied only to particular configurations: relatively low 



column densities (Richling et al. 2001; Richling 2003), Hubble 
flows (Kobayashi & Kamaya 2004), ID geometry (Ahn et al.), 
spherically symmetric configurations Dijkstra et al. 2005 ab); 
others are strongly tailored towards cosmological simulations 
(Cantalupo et al. 2005, Tasitsiomi 2006) and can deal with 
clumpy/inhomogeneous media (Spaans, 1996; Richling, 2003; 
Hansen & Oh 2006). In addition, the effect of dust absorption, 
which is one of the most important factor affecting Lya trans- 
mission, is treated only in some of the codes above. 

In any case, none of these studies has attempted to explain 
systematically the observed variety of Lya line profile mor- 
phologies. Furthermore, none of the above schemes has so far 
attempted a detailed modeling of individual galaxies, taking 
into account the available observational constraints, i.e. con- 
straints on the stellar populations, the ionised and neutral inter- 
stellar medium, on dust extinction and including their spatial 
distribution and kinematics. All this lack in the actual theoret- 
ical modeling needs to be filled up, in order to extract some 
information from the huge reservoir contained in the avail- 
able observational data. With these objectives in mind, and 
with the main aim of improving our understanding of Lya in 
both nearby and distant starburst galaxies, we have developed 
a general-purpose 3D Lya radiation transfer code applicable to 
arbitrary geometries and velocity fields. 

In the present paper we provide a description of the code 
and test its validity against known solutions and results from 
other codes reported in the literature. Exploring different ge- 
ometries, dust-free and dusty media, and different input spec- 
tra (e.g. line emission, or continuum + line), we examine 
the resulting line profiles and their dependence on various 
physical parameters. Our immediate goals are to obtain an 
overview over the possible Lya line profile morphologies, and 
to gain physical insight into the processes governing them. 
Applications to observed galaxies and other simulations will 
be presented later. 

The remainder of the paper is structured as follows. A de- 
scription of the radiation transfer code is given in Sect. 13 Tests 
of the code and results for simple geometrical configurations 
(slabs, infalling/expanding halos) are presented in Sect. [3] In 
Sect. |4] we comment on the formation of damped (Voigt) Lya 
profiles and related profiles. Spherically expanding dust-free 
or dusty shells are re-examined in Sect. [5] An overview of the 
predicted Lya line profile morphologies and qualitative com- 
parisons with observations is given in Sect. [6] Our main con- 
clusion are summarised in Sect.0 



2. Radiation transfer code 

A general 3D radiation transfer code MCLya allowing for ar- 
bitrary hydrogen density, ionisation & temperature structures, 
dust distributions, and velocity fields was developed using a 
Monte Carlo technique. The input files and the structure of the 
code have been designed for a future joint use with the 3D ra- 
diation transfer and photoionisation code CRASH of Maselli 
et al. (2003). We now summarise the main ingredients and as- 
sumptions made in this code. 
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2.1. Geometry 

The present version of the code assumes a 3D cartesian grid 
of n 3 cells. Typically we adopt n = 128. The relevant quanti- 
ties describing a 3D structure are the neutral hydrogen density 
distribution, the dust density distribution, the temperature dis- 
tribution, and the velocity field. These are prescribed by input 
files. 

2.2. Photon sources 

Lya and/or continuum photons are emitted from one or several 
point sources. Each source is described by : 

- its location, 

- the total number of emitted photons 

- optionally their emission direction, if not isotropic, 

- the source spectrum (typically monochromatic, a constant 
photon density per frequency or wavelength interval, a 
Gaussian, or combinations thereof). 

2.3. Physical processes 

To capture the essentials of radiation transfer in the UV in- 
cluding and around the Lya line we include three main physi- 
cal processes, dust absorption and scattering and the Lya line 
transfer, in the present version of our code. Given the princi- 
ples of Monte Carlo simulations, other processes can easily be 
included in the future, if desirable. 

2.3.1 . Lya line transfer 

In the whole section, we describe the Lya radiative transfer 
equations in a static medium. To adapt them to moving media, 
we just convert frequencies to local co-moving frequencies and 
convert them back to the external frame by a Lorentz transfor- 
mation, neglecting terms of order 0(v 2 /c 2 ). 

A Lya photon corresponds to the transition between the 
n = 2 and n = 1 levels of a hydrogen atom. This is the 
strongest H i transition, with an Einstein coefficient given by 
Ai\ = 6.265 x 10 8 s _1 . The scattering cross-section of a Lyo- 
photon as a function of frequency in the rest frame of the hy- 
drogen atom is: 



■fn 



r/4n 2 



m e c (v- v ) 2 + (r/47r) 2 ' 



(1) 



where fa = 0.4162 is the Lya oscillator strength, vo = 2.466 x 
10 15 Hz is the line center frequency, and r = An is the damping 
constant which measures the natural line width. 

The optical depth t v (s) of a photon with frequency v trav- 
eling a path of length s is determined by convolving the above 
cross-section with the velocity distribution characteristic of the 
absorbing gas, and is of the form: 



s*S poo 

r v (s) = I I m 

JO J-oo 



(V z )cr v dV z dl, 



(2) 



where V z denotes the velocity component along the photon's 
direction. Thermal motions of Hydrogen are described by a 



Maxwellian distribution of atoms velocities whose velocity dis- 
persion, V t h = (2Icb T/mu) 1 ^ 2 = 12. 85T^ 2 km s _1 , corresponds 
to the Doppler frequency width Av D = (V,h/c)vo. Here T4 is the 
gas temperature in units of 10 4 K. In certain cases an additional 
turbulent motion, characterised by V tur b, is taken into account 
in the Doppler parameter b given by 



b 



+ V; 



turb' 



(3) 



Let us now introduce some useful variables. First the fre- 
quency shift in Doppler units 



x 



■ 



V 



Av D b ' 



(4) 



where the second equation gives the relation between x and a 
macroscopic velocity component V measured along the photon 
propagation (i.e. parallel to the light path and in the same direc- 
tion). Second the Voigt parameter a = ^ = 4.7 x 1(T 4 T~ l/2 , 
or more generally a = 4.7 x 10 -4 (12.85 km s _1 /b) for non-zero 
turbulent velocity. Adopting this notation, it can be shown that: 

t x (s) = cr H (x)n H s= 1.041 x 10- 13 r 4 - 1/2 jV H (5) 

yn 

where tin is the neutral hydrogen density, and Nh the corre- 
sponding column density. The Hjerting function H(x,a) de- 
scribes the Voigt absorption profile, 



a f°° e-y 2 dy ( e^ 2 if 

H(x,a) = - ~ \ a Jf 

n J-00 (y -x) 2 + a 2 {-^ it 



\x\ < x c 
\x\ > x c 



(6) 



which is often approximated by a central resonant core and 
power-law "damping wings" for frequencies below/above a 
certain boundary frequency x c between core and wings. For 
a in the range of 10~ 2 to 10~ 6 , x c varies typically from 2.5 to 
4. To evaluate H in our code, we use the fit formulae given by 
Gray (1992). 

To characterise the depth of a static medium we will use To, 
the optical depth at line center: 



To(s) 



i.04i xio- 13 T: l/2 N H H{0 ' a) 



yn 



3.31 x 10~ 14 T~ ll2 N H 



for zero turbulent velocity, or more generally 
t (j) = 3.31 x 10~ 14 (12.85 km s-'/^jVtf 



(7) 
(8) 

(9) 



This monochromatic optical depth to has been used in most 
recent studies (e.g. Ahn et al. 2001, 2002; Zheng & Miralda- 
Escude 2002, Hansen & Oh 2006); however, it differs from 
the total Lya optical depth T^ euf used in the classical work of 
Neufeld (1990) by T^ euf = t V^- 

Once the absorption probability is given and before the 
Lya photon is re-emitted, its frequency and angular distribu- 
tion must be determined. If the atom is not perturbed by colli- 
sions during the time a Lya photon is absorbed and re-emitted, 
the frequencies before and after scattering are identical in the 
atom's rest-frame. On the other hand, when the atom undergo 
a collision, the electron is reshuffled on another energy level 
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and the frequencies before and after scattering are uncorrelated. 
Given the typically low densities in astrophysical media, we as- 
sume coherent scattering in our code. 

Concerning the angular redistribution, our code can model 
the case of isotropic as well as the more realistic dipolar re- 
distribution. In the case of isotropy, for simplicity and speed 
we use the angle averaged frequency redistribution function 
Rii-A(Xm, *out) from Hummer(1962). In practice we use pretab- 
ulated values of the cumulative frequency distribution function 
°f Rii-a f° r different input frequencies and temperatures. In 
all the static geometries presented in this paper, we use an 
isotropic angular redistribution. As test calculations confirm 
(cf. Fig^and also Zheng & Miralda-Escude 2002 and Hansen 
& Oh 2006) this is an excellent approximation since the Lya 
photons undergo a very large number of scatterings where any 
angular preference is smeared out. 

To avoid numerous core scatterings in static cases with a 
high column density, different acceleration methods have been 
developed in other radiation transfer codes (e.g. Ahn & Lee 
2002, Djikstra & Haiman 2005a). In our case, such an accel- 
eration is easily included in the redistribution functions by set- 
ting artificially Rn-A(xin, x ou t) = £ where e is close to zero for 
| jci n | < x c and |x out | < x c , which corresponds to setting the prob- 
ability to be re-emitted at |x out | < x c to zero, when photons 
are absorbed in the core (i.e. when |x; n | < x c ). In practice we 
have not used the acceleration method except for test cases, as 
it turns out that all cases shown here are tractable without it. 

The dipolar angular redistribution has been implemented 
without the use of redistribution functions, but microscopically, 
following the detailed descriptions of former codes (cf. Zheng 
& Miralda-Escude 2002 and Dijkstra & Haiman 2005a). The 
necessity to use this more physical redistribution is particularly 
important for expanding shells (see section[5j. 

2.3.2. Dust scattering and absorption 

During its travel in an astrophysical medium, the Lya photon 
will diffuse on H atoms, but it can also interact with dust: it can 
either be scattered, or absorbed. The dust cross-section <Td is 
composed of an absorption cross-section cr a , and a scattering 
cross-section <x v : 



cr„ + <r. 



(10) 



where cr a s — nd 1 Q a , s , with d the typical dust grain size which 
will affect Lya photons, and Q a s the absorption/scattering ef- 
ficiency. At UV wavelengths the two processes are equally 
likely, Q a * Q s * 1, so the dust albedo A = Q s /(Q a + Q s ) 
is around 0.5: half of the photons interacting with dust will be 
lost, and half will be re-emitted in the Lya line. 

We assume that the dust density rid is proportional to the 
neutral H density in each cell 



n d = n H X 



mn Md_ 
m d Mh' 



(ID 



where m d is the grain mass and m H the proton mass. The rel- 
evant quantity, r d given just below, is described by one free 



rrid = 3 X 10 17 g. The total (absorption + scattering) dust opti- 
cal depth seen by a Lya photon is then: 



Td = T a + T, 



-f 

JO 



o~d nd(s) ds. 



(12) 



The relation between the dust absorption optical depth at 
Lya wavelength t„ = (1 - A)tj and the color excess Eg-v is 
given by 



-B-V 



1.086- 



'1216 



R 



r a * (0.06... 0.1 1)t 



(13) 



where Ay and A\2\(, is the extinction in the V band and at 1216 
A respectively, R the total-to-selective extinction. The smaller 
numerical value corresponds to a Calzetti et al. (2000) attenua- 
tion law for starbursts, the larger to the Galactic extinction law 
from Seaton (1979). 

2.4. Monte Carlo radiation transfer 

For each photon source we emit photons one by one, and we 
follow each photon until escape from our simulation box or 
absorption by dust. Let us now describe one photon's travel. 

2.4.1 . Initial emission 

The emission of a photon is characterised by an emission fre- 
quency and direction. The frequency v (here in the "external", 
i.e. observer's frame) samples the source spectrum represent- 
ing usually Lya line emission and/or UV continuum photons. 
For media with constant temperature, v or more precisely the 
emission frequency shift from the line center, is conveniently 
expressed in Dopper units, i.e. x = (v - vo)/Avo (cf. Eq.@J. 

We assume that the source emission is isotropic (in the lo- 
cal co-moving frame, if the considered geometry is not static). 
Thus the emission direction, described by the two angles 8 and 
<f>, is randomly selected from 



6 = cos~\2fr - 1) 
(p = In 



(14) 
(15) 



parameter, the dust to gas ratio jf- assuming d = 10 b cm and 



A/n 



where £1,2 are random numbers l , < £1,2 < 1. The photon 
travels in this direction until it undergoes an interaction. In 
moving media, the photon frequency in the external frame is 
evaluated by a Lorentz transformation. 

2.4.2. Location of interaction 

The location of interaction is determined as follows. The op- 
tical depth, Tint, that the photon will travel is determined by 
sampling the interaction probability distribution P(r) = 1 - e~ T 
by setting 

7W = -Zn(l-£) (16) 

where £ is another random number. 

1 The random numbers generator used in the code is the ran func- 
tion from Numerical Recipies in Fortran 90 (Chapter B7, page 1 142) 
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We sum the optical depth r along the photon path s: 

T(S)=T X (S)+T d (s) (17) 

and we determine the length s corresponding to t(s) = t,„,. We 
calculate the coordinates corresponding to a travel of length s 
in the direction (6, 0) starting from the emission point. This is 
the location of interaction. Now, we have to compute if the Lya 
photon interacts with a dust grain or a hydrogen atom. 

2.4.3. Interaction with H or Dust ? 

The probability to be scattered by a hydrogen atom is given by: 
n H cr H (x) 



Pk(x) 



»ffO" H (x) + n d o- d ' 



(18) 



where cr H (x) = / K 



m "cAvd fl -* * s t ' le hydrogen cross section 
for a Lya photon of frequency x. We generate a random number 
< £ < 1 and compare it to P H : if £ < P H , the photon interacts 
with H, otherwise it is scattered or absorbed by dust. 



2.4.4. Scattering on H atoms 

When the photon is absorbed by an H atom, how will it be 
reemitted ? 

We first convert the frequency in the external (observers) 
frame, v eXt , to the comoving frequency of the fluid, v„„, with a 
Lorentz transformation 



Vim = V„,(l ), 



(19) 



where k, is the photon direction and V the macroscopic/bulk 
velocity of the H atoms. As already mentioned above, we as- 
sume partially coherent scattering and either isotropic or dipo- 
lar angular redistribution. After scattering, the new frequency 
is again converted back to the external frame. 

2.4.5. Dust scattering and absorption 

When the photon interacts with dust, we generate a random 
number < £ < 1 determining whether it is absorbed or scat- 
tered. In practice, if £, < A = Q s /(Q a + Qs) the photon is 
scattered by dust and simply reemitted coherently. Otherwise 
the photon is absorbed by dust and is considered lost for the 
present simulation. For the same reasons already discussed 
above (Sect. 12.3. 1> we currently assume that dust scatter- 
ing is isotropic, but other angular distributions (e.g. Henyey- 
Greenstein functions) can easily be implemented. 



. MC spectra 
Acceleration 

. Neufeld 
dipolar redist 




Fig. 1. Predicted emergent Lya profiles for monochromatic line 
radiation emitted in a dust-free slab of different optical depths 
(solid lines) compared with analytic solutions from Neufeld 
(1990, dashed). The dotted blue curve shows the line profile 
obtained using a frequency redistribution function, which skips 
a large number of resonant core scatterings. The adopted con- 
ditions of the medium are: T = 10 K (i.e. a = 1.5 10~ 2 ), and 
To = 10 4 , 10 s , 10 6 from top to bottom. The green long-dashed 
curve, obtained with a dipolar angular redistribution, overlaps 
perfectly the black solid line, obtained with the isotropic an- 
gular redistribution function, illustrating the fact that in static 
media, isotropy is a very good approximation. 



. MC spectrum 
. Neufeld 



Fig. 2. Same as Fig^foi" a slab of small optical depth (to = 
10 4 ) and high temperature T = 10 4 K(a = 4.7 10~ 4 ) illustrating 
the deviation from the analytic solution of Neufeld (1990) due 
to the more complete treatment of the line absorption profile 
(see text). 



2.4.6. Output 

The precedent scheme is repeated until the photon escapes the 
simulation volume, or undergoes a dust absorption. Then, we 
store all the information concerning this photon in a matrix and 
start with the next photon. This procedure is repeated for all 
photons and all emission sources. Finally all the desired re- 
sults, such as spatially integrated line profiles, monochromatic 
or integrated Lya images, surface brightness contours for any 



given line of sight and spatial resolution etc. are computed from 
this output matrix. For reasons of symmetry, and to maximise 
the numerical efficiency, all the spectra presented hereafter are 
integrated spectra over all directions, except if mentioned oth- 
erwise (cf FiglTTY 
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3. Validation of the code and examination of 
simple geometrical configurations 

To validate/test our code and to gain insight into basic prop- 
erties of the Lya radiation transfer we have computed a large 
number of simulations for geometrical setups discussed pre- 
viously in the literature : slabs, expanding/infalling halos, ex- 
panding shells that we will present in details hereafter, and 
disk-like configurations for which we have perfect agreement 
with Richling et al. (2003). 

We consider various input spectra (i.e. the intrinsic emis- 
sion line profile and possible continuum emission), especially 
the limiting cases of a pure monochromatic (i.e. line) radiation. 
Both cases with or without dust are considered. The case of a 
source emitting a pure continuum in a plane parallel slab with 
and without dust is discussed in Sect. [4] 

3.1. Homogeneous slab 

The best studied case is that of a plane parallel homogeneous 
slab, for which analytic solutions of the Lya transfer problem 
have been derived and discussed in depth by Neufeld (1990). 
For a given source position and input spectrum, the main phys- 
ical quantities determining the output spectrum are the temper- 
ature of the medium and its optical depth, whatever the angu- 
lar redistribution is. Note that in this section, all profiles have 
been obtained using isotropic angular redistribution, except if 
the contrary is specified. 

3.1 .1 . Monochromatic radiation, no dust 

As a first test we simulate a dust-free slab with a central 
plane source emitting Lye photons in the line center (i.e. vn = 
2.46 X 10 15 Hz or equivalently x = 0). We choose T = 10 K 
(i.e. a - 1.5 x 10~ 2 ), and vary the optical depth in the line cen- 
ter, to = 10 4 , 10 5 , 10 6 , to be in the validity range of Neufeld's 
analytic solution assuming a power-law absorption profile (Eq. 
|6j: a very optically thick slab, where a to > (a To) cr it (cf. be- 
low). In this case the emergent Lya line profile is given by (see 
Neufeld 1990, Eq. 2.24): 



/(t , x) - 



V6 



1 



24 aro CO sh[ ^[n ^ ^5A{\x i \|a^o)\ 



(20) 



As shown in Fig.^ our spectra are in perfect agreement with 
Neufeld's predictions 2 . The spectra are double peaked and 
symmetric around x — 0. The peak frequency x p reflects the 
physical properties of the neutral medium 



±0.88(flT ) 



1/3 



(21) 



As expected, the more optically thick the medium is, the 
more separated are the peaks. The width of the peaks becomes 
larger with larger tq. In Fig.n we a ls° show that computations 



2 As already mentioned above our definition of To, the monochro- 
matic line optical depth at x = (line center), differs from Neufeld's 
definition of T? euf which is the total, i.e. frequency integrated, Lya 



.without dust 
with dust : escaped photons 
.with dust : absorbed photons 




Fig. 3. Predicted emergent Lyo- profiles for monochromatic line 
radiation emitted in a slab of optical depth to = 10 6 , T = 10 K, 
and T a = 10~ 3 and 0. The black solid line shows the emergent 
profile without dust and the blue dotted line the one with dust. 
The red dashed line shows the (unobservable) spectrum of the 
photons which have been absorbed by dust. The effect of dust 
seems "grey" on the emergent line profile because Td ust is very 
weak. For greater Td us t, the inner part of the peaks is destroyed. 



optical depth. One has: r. 



■Neuf 



: T V^- 



using our acceleration method (cf. Sect. 12.3. It yield excellent 
agreement. 

For higher temperature, i.e. smaller a (see Fig. [5) our spec- 
trum differs from the analytic solution of Neufeld (1990), as 
also noticed by Ahn et al. (2001) and Zheng & Miralda-Escude 
(2002). This is due to the simplified assumption of a power-law 
line profile in the wings, which Neufeld (1990) assumes to be 
valid for the entire absorption line profile H (cf. Sect. 12. 3.11 . As 
a consequence our peaks are more separated and less symmet- 
ric than expected from Neufeld's calculation. This is due to the 
fact that Neufeld's approximation of H underestimates the ab- 
sorption probability in the core, such that photons escape more 
easily than in the real case and their mean escape frequency 
remains closer to the line center. 

For Neufeld's analytical solution, based on the assumption 
that H(x, a) » —Srr? as in Lorentz wings, to be valid, a mini- 
mum criterium is that the optical depth at the frequency x c cor- 
responding to the transition between core and wings is larger 
than 1 . Otherwise, most of the Lya photons will escape from 
the core where Neufeld's approximation is not valid. In prac- 
tice, since T cr i t = t(x c ) = 1 corresponds to (aTo)ait = V^Jy, 
the analytic solutions of Neufeld (1990) are valid only for suf- 
ficiently large optical depths, i.e. for a to > (a To) C rit- 

3.1.2. Monochromatic radiation and dust 

We now include dust in the slab and examine its effect on line 
profiles and on Lye photon destruction. 

As illustrated in Fig [3] the presence of dust absorption re- 
duces the total Lya intensity. For the case of monochromatic 
radiation with a very small amount of dust shown here, the 
shape of the emerging Lya profile remains basically unchanged 
by dust. However, with increasing Tdust> the inner parts of the 
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Fig. 4. Escape fraction f e for Lya line photons emitted 
monochromatically in a dusty slab as a function of (aro)^ 3 r a . 
The excellent agreement between our simulation (crosses), the 
computations of Hansen & Oh (2006) and the approximate an- 
alytical solution of Neufeld (1990) given by Eq.|23]is shown. 



emission peaks are destroyed. Although not an observable, we 
also plot the spectrum of the absorbed photons (red dotted line 
on Fig[3}. As expected, it is symmetric with respect to the line 
center: it presents two peaks, closer to the line center than those 
of the emerging profile. To understand these features, we have 
to take into account two competing effects, the large number 
of scatterings photons undergo in the line core favouring a po- 
tential interaction with dust, but the low probability to interact 
with dust in presence of the very strong HI absorption cross 
section close to the line core. Indeed, the probability to interact 
with H is around unity in the core since 



Ph(x = 0) = 



ct-h(O) 



°h(0)+ — j^crj 



l+c , T mM± 

4 M H 



-l 



(22) 



with c' = 3.3 10~ 6 Kr 1 / 2 for the dust parameters adopted here. 
Therefore none of the photons are absorbed by dust in the very 
center of the line. The peaks are located where the probability 
to interact with H is lower (for |;t| < x c ) but where the number 
of interactions (which ultimately increase the chances for dust 
absorption or scattering) is still high, i.e. on both sides of the 
line core. 

For an optically thick slab (ot > 10 3 ), Neufeld (1990) 
gave an analytic solution to determine the escape fraction, i.e. 
the ratio between the number of photons which escape the 
medium and the total number of emitted photons. This escape 
fraction f e depends on the combination {aTo) l ^T a , where T a is 
the optical depth of absorption from the center to the surface 
of the slab, i.e. t„ = t<//2 here since we consider a dust albedo 



A ~ 0.5. For a central source and in the limit (aro) 



1/3 



» T a 



Neufeld (1990, his Eq. 4.43) has derived an approximate ex- 
pression for f e , which, in our notation, is 



f e = 1 / cosh 



V3 



127 M'N 



1/2 



(23) 



where £ m 0.525 is a fitting parameter. In Fig.@J we compare 
our results with this analytical curve and with the results from 



Monte-Carlo simulations of Hansen & Oh (2006). Our results 
are in good agreement with both of them. 

3.1 .3. Continuous input spectrum 

The case of a source emitting a pure continuum in a plane par- 
allel slab with and without dust is discussed in Sect.@] 

3.2. Expanding/infalHng halos 

We now simulate spherical clouds ("halos") of uniform density, 
not only static, but also expanding or collapsing. An isotropic 
Lya source is either located at the center of the sphere, or uni- 
formly distributed over the whole volume. In this section, the 
angular redistribution considered is isotropic, after verification 
that the emergent profiles are very similar, even in the cases 
with motion. 



3.2.1. Monochromatic source 

We first chose the same physical conditions as Zheng & 
Miralda-Escude (2002) and performed a run with a colunm 
density measured from the center to the edge of the cloud 
Nu = 2. 10 20 cirT 2 , corresponding to a line-center optical depth 
To = 8.3 10 6 and a temperature of 2. 10 4 K. The velocity field 
is a Hubble-type flow, V(r) = +V m . dx r / 7? max , with V max = 200 
km s at the outer radius /? max of the halo, and the positive 
(negative) sign corresponding to the expanding (infalling) halo. 
Our results, shown in Fig.|5] are in good agreement with those 
of Zheng & Miralda-Escude (2002). For static halos, the same 
characteristic double peak profile as for the static slab (cf. Fig. 
[0 are obtained both for a uniform emissivity or for a central 
point source. As shown by Djikstra & Haiman (2005a) the po- 
sition of the peaks shows the same dependence as for the static 
slab (Eq. l2U . but is given by x p = + 1 .06(aro)^ 3 for oto > 10 3 . 

Lya line profiles from an expanding or infalling halo are 
perfectly symmetric to each other (compare the red short- 
dashed and blue long-dashed curves on Fig.|5}. Expanding ha- 
los present a red peak, whereas infalling halos have a blue one. 
Halos with a uniform emissivity (right panel in Fig. [5} show 1 ) 
broader lines, 2 ) emission extending on both sides, and 3) a sec- 
ondary peak on the blue (red) side for the expanding (infalling) 
halo. The last feature is not visible on the plots of Zheng & 
Miralda-Escude (2002), but their resolution may be too low. 
However, the results of Dijkstra & Haiman (2005a) show this 
secondary peak. We now briefly discuss the origin of these fea- 
tures. 

Why are single peaks formed in expanding/infalling media 
with a central point source emitting monochromatic radiation 
at the Lya line center? The reason is simple. The probability 
to escape the medium for a photon at line center is e~ T °, i.e. 
close to zero for both cases shown here. As an expanding halo 
contains atoms with velocities v(r) from to Vmax, all photons 
in the frequency range x = [0, V mlix /b ~ 11] will be seen in 
the line center by atoms of the corresponding velocity, and are 
thus "blocked". Therefore the only possibility to escape is to be 
shifted to the red side. The symmetry of the double peak profile 
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Fig. 5. Predicted Lya profiles of an isotropic and monochromatic source in uniform and isothermal spherical halos with a center 
to edge column density of Nh - 2. 10 2() cirT 2 and T - 20000 K (i.e. b = 18.24 km s _1 ). Black solid lines show the Lya profiles 
for a static halo, red short-dashed lines the expanding case with V max = 200 km s -1 , and blue long-dashed the infall case (same 
V max ). Left: Central point source. Right: Uniform emissivity. On the right panel, the "infalling" spectrum has been decomposed 
in two components: the emergent spectrum from an "inner part" (the inner 50 % of the volume), resembling that from a halo 
with a central source, and the one from the "external part" akin to the spectrum from an expanding shell where photons do not 
undergo backscatterings (cf. Fig. I12i. 




Fig. 6. Predicted Lya profiles of an isotropic and monochromatic source in uniform and isothermal spherically expanding halos 
with varying column density (Nh = 2. 10 18 to 7. 10 20 cirT 2 as indicated in the Figure), and other parameters fixed to T — 20000 
K (i.e. b = 18.24 km s -1 ) and V max = 200 km s _1 . Left: Central point source. Right: Uniform emissivity. 



of the static case is "broken" in this way and transformed to a 
red peak for an expanding halo. 

The increased line width in case of a spatially uniform 
emissivity (1) and the presence of photons on both the red and 
blue sides (2) is due to the fact that the intrinsic line emis- 
sion (assumed at x — in the atom's frame) spreads already 
over all frequencies from x - -V m:lx /b to +V max /£> in the ex- 
ternal (observer's) restframe. Radiation transfer effects further 
redistribute the photons in the wings. In fact, for the expand- 
ing halo, photons emerging with x > V max /£> (i.e. very blue 
ones) correspond to photons emitted close to the halo edge on 
the approaching side, which have been further redshifted by 
diffusion away from the line center (at x — HV max /£> in the ob- 
servers frame). This naturally also produces the local minimum 
observed at this frequency, which separates the secondary peak 



from the main one (point 3 above). This is easily verified by 
plotting e.g. the contributions from photons emitted in the "ex- 
ternal" parts (i.e. close to the edge) and from "internal"regions, 
as shown in the right panel of Fig.|5] We chose as inner/outer 
limit the radius which corresponds to an inner sphere and an 
outer shell of same volume. Therefore the location of this min- 
imum is a measure of the external velocity of the halo, as al- 
ready noticed by Dijkstra & Haiman (2005a), at least for ve- 
locities V max ^ 360 km s _I beyond which the two peaks start to 
overlap. 

With Figs.|S]and0 we investigate how the emergent lines 
depend on column density and velocity gradient. When in- 
creases, the emission peak is shifted away from line center and 
broadens out for both cases (central point source and uniform 
emissivity). Indeed, the optical depth at line center increases 
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Fig. 7. Predicted Lya profiles from an expanding halo with a 
central monochromatic source and varying velocity gradients 
(y max = 20, 200, and 2000 km s" 1 ) for T = 20000 K (i.e. 
b = 18.24 km s _1 ) and N H = 2. 10 2() cnT 2 . 

with jVh, so Lya photons have to diffuse far in the wings to es- 
cape. From the right panel of Fig.[6]we note that the location 
of the flux minimum between the two peaks remains constant 
quite independently of A^h- The behaviour for different veloc- 
ity gradients, first discussed by Wehrse & Peraiah (1979),is the 
following (cf. Fig. [7}: from V max = 20 to 200 km s _1 one goes 
from a nearly static case (i.e. double peaks) to a broad asym- 
metric line, whose peak position is progressively displaced red- 
wards. Above a certain value of the velocity gradient (or equiv- 
alently V max X the peak position moves back closer to line cen- 
ter as the medium becomes more "transparent" there, given 
the finite and constant column density. In any case, the ex- 
tent of the red wing continues to increase with increasing V max . 
Qualitatively the simulations shown in Figs. [6] and resem- 
ble those of a cosmological Hubble flow (cf. Loeb & Rybicky 
1999) although modified by a finite outer boundary. 

For test purposes we have also compared the surface bright- 
ness profiles with the results of Zheng & Miralda-Escude 
(2002), and find a good agreement. 

3.2.2. Continuous input spectrum plus a line and dust 

Considering a source spectrum with a pure continuum and 
varying amounts of dust yields the following for the same con- 
ditions considered before(see Fig. [8}. Due to the scattering na- 
ture of Lya, the continuum photons are removed from line cen- 
ter and redistributed in the wings leading to a P-Cygni type pro- 
file with a red emission and blue absorption in the expanding 
medium, and an inverted P-Cygni profile in the infalling case. 
In other words, a flat source spectrum can result in a complex 
line profile. 

The peak is located at the same frequency as in the 
monochromatic case with identical physical conditions, and 
the absorption feature is around x = V max /b. Again its loca- 
tion indicates the external velocity of the medium. The emis- 
sion peak is broader for uniform emissivity than for a central 
point source, as already discussed in the monochromatic case, 



whereas the absorption feature behaves inversely. The latter 
behaviour is due to a partial "replenishment" of the absorp- 
tion feature by the "secondary peak" discussed above for the 
monochromatic radiation. Also shown in Fig.0is how the pres- 
ence of dust changes the observed line profile: as dust absorp- 
tion is efficient in the wings close to (but not at) line center (cf. 
Sect. 13. 1 .21 it easily suppresses the emission peak and broad- 
ens the absorption feature leading to an asymmetric absorption 
profile with a more strongly damped red (blue) wing for an in- 
falling (expanding) halo. 

Far from the line center, the escape fraction, i.e. the ratio 
between the number of photons which escape the medium in 
one frequency bin and the number of photons emitted in this 
frequency bin, follows the expected exponential law f e = e~ T °, 
where t„ is the dust absorption optical depth. Indeed, when the 
influence of hydrogen scattering becomes negligible, photons 
are be absorbed by dust with the probability P = 1 - <?~ T ". 

To illustrate the variety of line profiles predicted for vari- 
ous intrinsic line strengths, we show in Fig.|9]the case of an ex- 
panding halo with uniform emissivity. As seen from this figure, 
a family of line profiles with intermediate cases between a pure 
continuum (cf. Fig. [8} and pure line emission (Fig. |5} is ob- 
tained with increasing intrinsic Lya equivalent width W(Lya). 
Note that for a sufficiently large W{Lya) the secondary peak be- 
comes again visible. Also, remember that for the case of oppo- 
site movement, i.e. infalling, the predicted spectrum is identical 
except for a change of red and blue frequencies (i.e. change x 
to -x). 

Overall, as we have seen, there are several degeneracies 
which make it difficult if not impossible to determine physi- 
cal parameters such as A^h and V max for cases of expanding or 
infalling halos. This is in particular complicated by the lack of 
a priori knowledge of the spatial distribution of the emissiv- 
ity and of the precise velocity field. For a sufficiently extended 
Lya emissivity it may be possible to determine the outer veloc- 
ity, V max , from a local minimum of the flux in the blue part of 
the line profile, as already pointed out by Dijkstra & Haiman 
(2005a). However, detecting this feature would require a fairly 
high signal to noise. Furthermore this blue side of the Lya pro- 
file may be altered by intervening IGM absorption components 
at slightly lower redshifts. 

4. On the formation of damped Lya absorption 
profiles 

As damped Lya profiles are frequently observed and pure Voigt 
absorption line profiles often used to fit components of Lya it 
is of interest to examine under which conditions actually pure 
(damped or non-damped) Lya absorption line profiles are ex- 
pected. 

To illustrate the point, and the non-trivial problem of Lya 
radiation transfer effects we show in Fig.[TO]the predicted spec- 
trum around Lya for an infinite slab illuminated uniformly in 
the central plane by a pure continuum source. Indeed from the 
discussion of the previous cases it is not surprising to find a 
double peaked emission profile with a deep central absorption, 
as shown for the case without dust. Adding already quite small 
amounts of dust allows one to destroy the emission peaks and 
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Fig. 8. Left: Line profiles for an expanding halo with a central monochromatic source showing the variety of profiles from P- 
Cygni to somewhat asymmetric damped profiles. Right: Infalling halo with a uniformly distributed emissivity showing profiles 
from inverted P-Cygni to less broad absorption. In both panels the adopted conditions are the following: Nh = 2. 10 20 cm -2 , 
T = 20000 K, and V max = 200 km s _1 . The line profiles are shown for dust amounts corresponding to: (Mj/Mh, T a ) = (0., 0.), 
(8. 10~ 4 , 0.1), (4. 10~ 3 , 0.5), (8. 10~ 3 , 1.). 
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Fig. 9. Predicted Lya- line profiles from an expanding halo 
with uniformly distributed sources emitting a flat continuum 
plus a monochromatic line (at Lya line center) with differ- 
ent equivalent widths. The physical conditions in the halo are: 
Nh = 2. 10 20 cnr 2 , T = 2 x 10 4 K, V max = 200 km s _1 , and 
no dust. The Lya line equivalent width, given in units of the 
Doppler width Avd, increases from to oo from the top-left 
corner to the bottom-right. 

to create Voigt-like absorption profiles, as also shown on Fig. 
ITOl All these profiles can be seen as the sum of two compo- 
nents: first photons which have not undergone any interaction 
in the slab leading to a Voigt profile, added to the usual double- 
peaked profile arising from photons which scattered in the slab. 
The latter one becomes less important when the amount of dust 
increases (for a fixed Nh) and/or when A^h increases (for a 
fixed dust-to-gas ratio). Hence the resulting profile approaches 
a Voigt profile. However, one should notice that the column 
density derived from Voigt-profile fitting to these predicted pro- 
files can be several times smaller than the "true" column den- 
sity in the simulation. For example, for the most dusty simula- 



tion (green long-dashed curve) from Fig. ^| the Voigt profile 
yields a good fit (in the sense of a reduced^ 2 ~ 1.02), and the 
best fit column density is a factor ~ 3.1 lower than the "true" 
input value. 

The reason for this apparently "strange" behaviour is sim- 
ply due to the fact that the predicted spectrum is computed 
here from all the emergent photons integrated over the whole 
object and over all emergent directions, i.e. it corresponds to 
an integrated spectrum of a symmetric object (e.g. the sphere 
with a central source) or of a sufficiently extended "screen" 
between us and the source. In this case we include all photons, 
some of which have undergone a complex scattering history be- 
fore emerging toward the observer. Therefore the double-peak, 
characteristic of static media is unavoidable without dust and 
the only way to "destroy" them is by adding dust. 

Considering different geometries such as a finite slab, or a 
sphere of the same center-to-edge optical depth and same tem- 
perature yields exactly the same line profiles as the infinite slab. 

The "implicit" geometrical assumptions made above for the 
construction of an integrated spectrum are unlikely to be appli- 
cable to the typical damped Lya systems (DLAs) observed in 
spectra of distant galaxies and quasars. These cases are better 
idealised by a small (in angular size) intervening cloud between 
the background source and the observer, which diffuses basi- 
cally all photons out of the line of sight. In other words pho- 
tons far from the line center will simply travel freely through 
the cloud, but when the photon frequency x is close to the line 
center, the probability to cross the medium without interac- 
tion, exp(-ToH(x, a)), is considerably reduced. Since photons 
which interact with the medium have basically no chance to 
be re-emitted along the same line of sight (compare the ob- 
server solid angle with An), they will be lost for the observer. 
Therefore the observed profile is simply that of a pure ab- 
sorption line described by the Voigt absorption line profile, 
namely I(x) = Iq exp(-To H(x, a)), reflecting the properties of 
the medium, i.e. the total column density and the temperature. 
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Fig. 10. Emergent Lya profiles for varying amounts of dust pre- 
dicted for a slab with an isotropic source in the center emitting 
a flat continuum. The conditions of the medium are: to = 10 4 , 
T = 10 K. In this static configuration, the two angular redis- 
tributions are equivalent; spectra have been obtained with the 
isotropic redistribution. The line profiles are shown for absorp- 
tion dust optical depths of r a = 0. (no dust), 0.1, 0.5, 1. from 
top to bottom. Note the double-peak profile plus central absorp- 
tion characteristic of a static medium for no dust. Already small 
amounts of dust are able to suppress the peaks, thereby leading 
to Voigt-like profiles. However, in such geometrical configura- 
tions, classical Voigt-profile fitting will underestimate the true 
HI column density. See text for more explanations. 



We simulated this configuration, and the resulting spec- 
trum is presented on Fig ^2 by the black curve. We checked 
again that a spheric cloud or a slab perpendicular to the line 
of sight, with the same optical depth and temperature, lead 
exactly to the same observed spectrum. For reasons of calcu- 
lation time, we considered a relatively small column density, 
N n - 2x 10 15 cra~ 2 . The temperature is set to T = 100K. As 
expected, the simulation is perfectly fitted by a theoretical Voigt 
profile (red dotted curve). Also plotted on the same graph is the 
reflected spectrum (blue curve), composed of photons which 
escaped the medium by the side they entered. 

Photons which are not "reflected", i.e. backscattered to the 
source will diffuse in the medium and finally escape after a 
large number of scatterings forming the "diffuse spectrum". 
For obvious reasons this spectrum, shown by the green curve, 
presents the same shape as the emergent spectrum from a slab 
with an isotropic source in the center, i.e. two symmetrical 
peaks. In dust-free cases with a flat incident continuum source, 
the transmitted spectrum will be the opposite of the reflected + 
diffuse spectrum, as all Lya photons are conserved. 

In short, to form pure symmetric Lya absorption line pro- 
files from a flat continuum source requires specific geometrical 
configurations allowing the photons to diffuse out of the line 
of sight. An alternative way to achieve such profiles, e.g. in 
an integrated spectrum of an embedded source, is by invoking 
the presence of dust which destroys the double emission peaks 
otherwise present. However, in such a case the apparent column 



Fig. 11. Predicted Lya profiles of the direct transmitted (black 
solid line), the reflected (blue short-dashed), and the diffuse 
(green dotted) radiation from a static slab illuminated perpen- 
dicularly by a distant background source with a flat continuum. 
The properties of the intervening slab are: N-r = 2. 10 15 cm -2 , 
T = 100 K (i.e. b — 1.29 km s _1 ). The directly, i.e. perpendic- 
ularly, transmitted flux yields the expected Voigt profile (red 
dotted line). The diffuse (scattered) light shows again the dou- 
ble peak profile as the slab with a central source. 

density derived from simple Voigt profile fitting underestimates 
the true value of Nu due to radiation transfer effects. 

5. Lya transfer in expanding dust-free and dusty 
shells revisited 

There are numerous indications, both theoretically and obser- 
vationally, for the presence of expanding shells and bubbles in 
starbursts. It is therefore important to simulate such geometri- 
cal configurations to examine both qualitatively and quantita- 
tively the diversity of Lya line profiles and to gain basic insight 
into the physical processes shaping them. 

Our model of an expanding homogeneous shell is described 
by the following parameters: an inner and outer radius R m i n and 
Rmax respectively, a uniform radial expansion velocity V exp , the 
radial colunm density Nu, and a constant temperature T given 
by the Doppler b parameter (Eq. [3}- The interior of the shell 
is assumed to be empty, the isotropic Lya/continuum emitting 
source located at the center. In contrast to the geometrical con- 
figurations discussed earlier, emergent profiles from expanding 
shells are sensitive to the angular redistribution. Therefore the 
dipolar redistribution is taken into account to treat the Lya ra- 
diation transfer consistently in this case, and all spectra shown 
in this section have been obtained using the dipolar redistri- 
bution in the code. The main parameters determining the Lya 
photon escape and hence the line profile are Vexp, Mi, b, and 
the dust-to-gas ratio M^/Mu, as shown below. The precise val- 
ues of the shell radii and thickness, setting the geometrical size 
and curvature, are secondary. 

Below we shall examine the following cases: The academic 
cases of a source with monochromatic emission and shells both 
without and with dust. These first two cases are essential to un- 



12 



Verhamme, Schaerer & Maselli: 3D Lya radiation transfer 




Fig. 12. Left: Emergent Lya profile from an expanding shell with central monochromatic source. The different shapes can be 
described with the number of back-scatterings photons undergo : bumps la and 2 are built-up with photons which didn't undergo 
any backscattering, the highest peak located at x = -2V exp /b (feature lb) is composed of photons which undergo exactly one 
backscattering, and the red tail Ic is made of photons which undergo two or more backscatterings. See text for more details. 
Right: Scheme of an expanding shell showing the origin of photons of different groups. The observer is on the left side, at 
infinity. 



derstand more realistic simulations allowing for arbitrary input 
spectra (including continuum and/or Lya line emission), again 
dust-free or with dust, discussed subsequently. 

5.1. Monochromatic emission and dust-free shells 

5.1 .1 . Basic line profile formation mechanism 

First we discuss the emergent Lya spectrum from a dust-free 
expanding shell with the following parameters: V exp = 300 km 
s -1 , N H = 2 x 10 20 cm- 2 , b = 40 km s _1 (i.e. r = 3.8 x 10 6 ), 
^min = 0.9/? max , and /? max = 1.17 x 10 19 cm. For monochro- 
matic Lya photons emitted at line center the resulting line pro- 
file is shown on the left panel of Fig ^] (solid black line). 
Qualitatively this line profile and others shown below exhibit 
the following characteristic features, marked on the figure and 
illustrated on the right panel: 

labc) an extended redshifted emission with one or two 
"bumps" (la and lb) and a red tail (lc) - all at x < (i.e. 
red side of Lya), 

2) a smaller blue bump (at x > 0), and 

3) an narrow emission peak at the line center (x = 0). 

Although we considered the same physical conditions as Ahn 
et al. (2003, their Fig. 2), features 2) and 3) are not apparent in 
their simulations (cf. also Ahn 2004). However, these features 
are also found in the simulations of Hansen & Oh (2006) and 
the origin of all of them is well understood, as we shall now 
discuss. To do so it is instructive to group the emergent pho- 
tons and to distinguish the emergent line profiles according to 
the number backscatterings they have undergone (see left panel 
on Fig. ll2t . A photon is said to "backscatter" when it travels 
across the empty interior before re-entering the shell at a dif- 
ferent location 3 (see right panel of Fig. 1121 . Note, that any such 

3 This definition is equivalent to the one used by Ahn et al. (2003). 



travel is counted as a backscattering, irrespective of its precise 
direction/length. In particular this does not necessarily imply a 
hemisphere change for the photon. 

Features la and 2: photons with zero backscattering. In 

this simulation all photons are emitted at line center (x = 0) at 
the center of the shell. Once a photon reaches the shell for the 
first time it is seen redshifted to -xcMF(V exp ) = -7.5 by the H 
atoms (in the comoving frame, CMF). A fraction of the photons 
will diffuse progressively through the shell towards the exterior 
and escape without backscattering (solid lines on the left panel 
of Fig.ll2i. Their spectrum (marked as backscattering) gives 
rise to an asymmetric double peak with a small blue component 
centered at X2 ~ 10 (feature 2 above) from photons escaping 
the blue wing of Lya in the (blueshifted) shell approaching the 

observer, and a somewhat redshifted stronger peak at x\ a 5 

(bump la above) corresponding to the photons escaping the 
red Lya wing in the blueshifted shell. Qualitatively this part 
of the spectrum is equivalent to the spectrum of a slab with 
a constant receding macroscopic velocity with respect to the 
emitting source (see Neufeld 1990, Fig. 6). 

Feature 3: direct escape. For sufficiently small column 
densities and/or large expansion velocities, a non-zero frac- 
tion of photons traverses directly the shell without interacting. 
This case appears with the probability e~ T ™p(°) ? where r exp (x) = 
t(x - x(V exp )) is the reduced Lya optical depth seen by a pho- 
ton with observers frequency x due to Doppler shift of V exp . For 
the case discussed here T exp (0) = t(— 7.5) = 6.09, so that 0.25 
% of the photons will escape without interacting. These pho- 
tons give raise to feature 3 in the bin at X3 =0 labeled in Fig. 
IT2l The importance of these direct photons increases of course 
with increasing V exp and decreasing column density, as seen in 
Figs. ll4l andll6l Before comparison with observed spectra this 
peaked flux contribution must obviously be convolved with the 
instrumental resolution. 



Verhamme, Schaerer & Maselli: 3D Lyo- radiation transfer 



13 



Features lb and lc: photons undergoing one or more 
backscattering. Let us now examine the situation after one 
scattering in the expanding shell. After the first scattering the 
frequency distribution of the photons in the external frame is 
shown by the black solid line on Fig. ^] As the photons are 
initially seen by the atoms at the frequency -x(V exp ), i.e. in the 
wings, they are basically re-emitted at the same frequency in 
the atoms frame. Depending on their emission direction this 
leads to a range of frequencies in the observer frame reaching 

from x 2x(V exp ) to 0, with more photons re-emitted in the 

absorption direction (around x 2jc(V exp )) or in the opposite 

direction (x ~ 0) because of the dipolar angular redistribution 
(the frequency distribution for isotropic redistribution would be 
nearly a square profile over x ~ [-2jc( V exp ), 0]). Which of these 
photons are now able to escape the medium after just one scat- 
tering is represented by the blue dotted histogram. 

Overall, depending on their frequency being in one of the 3 
spectral regions indicated in Fig.^J the fate of the photons is 
as follows. 

- Frequency range 1: Photons with a frequency x Q \, s ~ can- 
not escape the medium: they are emitted outward in the ra- 
dial direction with a frequency too close to the line center 
(in the atoms frame). Their escape probability is negligible. 

- Range 2: Around x b s 7.5, although their frequency is 

far from line center, no photons escape the medium for ge- 
ometrical reasons: their emission direction is perpendicular 
to the radial direction increasing thus strongly the geomet- 
rical path and the optical depth. 

- Range 3: Most of the photons escape with a frequency 
Xobs ^ -9 since their frequency is very far from line center, 
and their emission direction is convenient: they cross the in- 
ner part of the shell (backscattering), and when they arrive 
on the other side, the combination of their frequency and 
direction with the local macroscopic motion favours their 
escape. The frequency distribution function after one scat- 
tering shows that the number of photons re-emitted with a 
frequency x < -2jt(V exp ) decreases very rapidly. Therefore 
the photons escaping after just one scattering already show 
a peak close to the frequency corresponding to twice the 
expansion velocity. So a peak of escaping photons centered 
at the frequency x = -2jc(V exp ) appears. 

Photons undergoing further scatterings will be absorbed 
again, and the escape of those re-emitted around -2x(V exp ) 
will be favoured again for same reasons. This explains why 
the most prominent feature (lb) in the red part is located at 

Xib 2x(V exp ), measuring therefore twice the shell velocity. 

Photons undergoing progressively more scatterings will show 
a broadening frequency distribution compared to that after one 
scattering. The broadening of its red wing is responsible for 
the last feature (lc) made of photons escaping after 2 or more 
backscatterings. 

5.1.2. Dependence on shell parameters 

Another way to understand the different features of the emer- 
gent Lya line profile of an expanding shell is by varying the 
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Fig. 13. Frequency distribution of the photons in the expanding 
shell after the first scattering. The black solid curve contains all 
photons, the blue dotted one represents the histogram of pho- 
tons which escaped after only one scattering, bf They form a 

bump around x 2x(V exp ) which explains the appearance of 

feature lb. See the text for more details. 
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Fig. 14. Dependence of the predicted Lya profiles from an ex- 
panding shell with a central monochromatic source on the ex- 
pansion velocity. 

parameters. Let us examine how the spectrum depends on the 
expansion velocity V exp , on the thermal and turbulent veloci- 
ties intervening in the Doppler parameter b, and on the colunm 
density A^h- 

For increasing expansion velocities, and starting from the 
static case with a symmetric double peak profile (cf. above 
and Eq. B18 of Dijkstra & Haiman 2005a for a static sphere), 
the imbalance between the red part and the blue part grows 
(see Fig. I14l> : progressively more photons escape from the red 
part of the line because atoms see them already redshifted (red 
dashed curve) at the first interaction. The probability to be re- 
emitted in the line core is then smaller than the one to "remain" 
in the wing. Hence the growing asymmetry between red and 
blue. Note also the appearance of excess flux at line center for 
the curve with the highest V exp in Fig.^] The appearance and 
strength of this feature (xj, = 0, feature 3) is consistent with the 
increasing direct escape probability. 
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Fig. 15. Dependence of the predicted Lya profiles from an expanding shell with a central monochromatic source on the Doppler 
parameter b. As the characteristic unit of Lya radiative transfer problem, x, depends on the Doppler parameter b(cf. Eq. 0J, a 
variation of b will lead to a change of x too. To compare cases with different b, we show the usual plot in x units but also a plot 
in velocity. Left: Plot in Dopper units x, Right: Plot in velocity along the line of sight v D b s . 



For a given b, the red part is a single peak, for small enough 
values of V eX p (e.g. V exp < 200 km s _1 for b — 40 km s _1 ). In 
fact, the two contributions of photons undergoing zero and one 
backscattering are then too close to be distinguished, and the 
two corresponding red peaks (features la and lb) are mixed. 
The same trend is found for large values of b, as shown in the 
left panel of Fig. ^] More precisely, predicted Lya profiles 
from an expanding shell present only one red peak for suffi- 
ciently small V exp /b. With increasing V exp or decreasing b, the 
blue part (feature 2) becomes very faint, and almost all photons 
escape in the red part, presenting two well separated peaks. In 
all cases, and quite independently of b (cf. Fig. [21 left), the 
second and most prominent peak (feature lb) traces twice the 
expansion velocity as already discussed above. It is essentially 
composed of photons which undergo one backscattering, while 
the first peak (feature la) is made of photons undergoing no 
backscattering. Note, however, that in contrast to the appear- 
ance in this plot, the position of main peak is independent of b 
in observed spectra, as illustrated in the right panel of Fig. 1151 
This is simply due to the definition of x, which depends on b. 

Varying the column density (leaving the other parameters 
unchanged) has the following effects (see Fig 1 1 61 . First, the 
flux excess at X3 = 0, i.e. the fraction of photons which can es- 
cape the shell without interacting, decreases exponentially with 
increasing Nh- Second, the relative importance of the two red 
peaks (la, lb) changes: with increasing column density the first 
peak (feature la) decreases with respect to the second one (lb), 
whereas the red wing (lc) is enhanced, since the importance of 
backscattering increases. For sufficiently large Nh the first peak 
(la) disappears completely, whereas the red tail (lc) becomes 
as important as the mean peak (lb). 

In short, the imbalance between blue and red emission in- 
creases with increasing V exp /b. The separation between the 
multiple peaks formed on the red side of the Lya profile be- 
comes progressively less clear (i.e. the peaks merge together) 
for lower expansion velocities, and/or higher temperature or 
turbulent velocities. For sufficiently large H 1 column densities 
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Fig. 16. Dependence of the predicted Lya profiles from an ex- 
panding shell with a central monochromatic source on the col- 
umn density. 

(Nh ~ 10 20 cirT 2 ) the main red emission peak measures quite 
well 2 V exp . 

5.2. Monochromatic emission and dusty shells 

In Fig. El we present the influence of dust on the emergent 
spectrum from the expanding shell studied above (Fig I12> . 
When t„ increases, the relative height of the two red peaks 
(features la and lb) is reversed: in dust-free media the promi- 
nent peak is lb, but when dust is present la becomes as high 
as lb. In the most dusty cases, corresponding to a destruction 
of ~ 93% of Lya photons, one notes a loss of photons from 
the the blue bump (feature 2 on Fig. I12i and from the red tail 
(feature lc on Fig ll2b . This is easily understood as these fea- 
tures are composed of photons undergoing a very large number 
of scattering, which increases their chance to be absorbed by 
dust. Qualitatively our results show the same behaviour as the 
outflowing shell with holes or clumps modeled by Hansen & 
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Fig. 17. Predicted emergent Lya profiles from an expanding 
shell with a monochromatic central source with a column den- 
sity Nh = 2 x 10 20 crrT 2 , an expansion velocity V exp = 300 
km s _1 a Doppler parameter b — 40 km s , and different dust 
absorption optical depths T a =0., 0.1, 0.5, and 1.0. The area 
under each curve is normalised to unity, i.e. corrected for pho- 
ton destruction by dust, to highlight the relatively minor line 
profile changes. The blue bump (2) and the red tail (lc) are 
more destroyed by dust than other features. Overall the profile 
is globally preserved, but somewhat "sharpened". 

Oh (2006): strongly redshifted photons are suppressed by dust, 
whereas the spectral peaks are still visible. The result is a some- 
what "sharpened" line profile. Although this overall "sharpen- 
ing" trend is also found in the earlier simulations of Ahn (2004) 
our results differ quite strongly from theirs, as already men- 
tioned above. 

We now quantify the Lya photon destruction by dust and 
its dependence on the shell and dust parameters. For illustra- 
tion we have chosen similar conditions as those discussed by 
Ahn (2004), namely a shell with an inner and outer radius 
^min = 0.9/? max , an expansion velocity V exp = 200 km s , a H i 
column density between A^h = 10 19 and 10 21 crrT 2 , a Doppler 
parameter of b — 20, 40, and 80 km s _1 , and a monochromatic 
central point source. In contrast to Ahn (2004) we assume that 
there is no dust inside the bubble, i.e. at R < R m i n , as most of 
the dust is probably destroyed there. Furthermore, the presence 
of dust inside the shell is not compatible with inferences from 
the empirical Calzetti attenuation law (see e.g. Gordon, Calzetti 
& Witt 2003). The predicted Lya escape fraction as a function 
of the dust absorption optical depth r fl measuring, for a given 
AT H column density, different dust-to-gas amounts, is shown in 
Fig.EU 

As expected, the main dependence of the escape fraction 
f e is on A^h and x a . Due to the multiple resonant scattering of 
Lya photons on hydrogen and the concomitant increase of the 
photon path length, f e is considerably smaller than the simple 
dust absorption probability exp(— r a ), and the escape fraction 
decreases strongly with increasing A^h- 

In an expanding shell the Lya photon destruction by dust 
depends also to some extent on the gas temperature (or on the 
Doppler parameter b), although in a somewhat "subtle" way as 
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Fig. 18. Escape fraction of Lya photons from an expanding 
shell with V e xp=200km s~'for different column densities (A^h = 
10 19 , 10 20 , 10 21 crrT 2 ) and for different Doppler parameters 
(b =20, 40 and 80 km s _1 ). 

can be seen from Fig. \TQ Indeed, for large column densities 
(here A^h > 10 20 ) f e is found to depend little on b, for interme- 
diate values of A^h ~ 10 20 the escape fraction varies in a non- 
monotonous way with b, and for lower A^h values f e decreases 
with increasing b. This latter behaviour is opposite to the one 
in a static medium, where the escape fraction increases with b 
(cf. Fig.|4]or Eq.[23]with a oc b~ x and to °c b^ 1 ). The reason for 
this inverted dependence is basically due to the fact that in an 
expanding shell the quantity V exp /b determines the frequency 
at which the initially emitted photons are seen by the atoms in 
the receding shell. For increasing b, V exp /b decreases, so that 
photons reaching the shell are seen with a frequency closer to 
line center. Hence they will diffuse more and will have a higher 
probability to be absorbed by dust. On the contrary, when b de- 
creases, photons are seen in the wings, their chance to interact 
with dust decreases and so does f e . When they are sufficiently 
far from line center (V exp /b » x c , where x c is some critical 
frequency), the shell becomes transparent in Lya and the es- 
cape fraction approaches the minimum value given by the dust 
absorption probability exp~ T " (red upper curve on Fig.1181. 

5.3. Dust-free and dusty shells with arbitrary source 
spectra 

Now we present the emergent spectrum from an expanding 
shell when the input spectrum is a flat continuum and with 
different amounts of dust. In Fig.^5] the black solid line rep- 
resents the emergent spectrum from a dust-free shell. It is a 
P-Cygni profile, quite similar to the expanding halo with a cen- 
tral point source in Fig. [8] Again, Lya radiation transfer leads 
to the appearance of an emission peak and an absorption fea- 
ture in the emergent profile which did not exist in the input 
spectrum. The first remark is that this spectrum is less complex 
than in the monochromatic case: there is only one red peak, 
and no blue bump, due to radiation transfer of photons at all 

frequencies. The emission peak is located at x 2x(V exp ), as 

the highest peak of the monochromatic spectrum, and the ab- 
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Fig. 19. Predicted Lya profiles for an expanding shell for a cen- 
tral source emitting a flat continuum and with different amounts 
of dust. The adopted conditions are the following: Nh = 2. 10 20 
cm" 2 , T = 10000 K, b = 40 km s _1 , and V exp = 300 km s" 1 . 
The line profiles are shown for dust amounts corresponding to 
(M d /M H , t b ) = (0., 0.), (8. 10~ 4 , 0.1), (4. 10~ 3 , 0.5), (8. 10~ 3 , 
1.). 
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Fig. 20. Predicted Lya profiles from an expanding shell with a 
central source emitting a flat continuum plus a monochromatic 
line (at Lya line center) with different equivalent widths. The 
physical conditions in the shell are: a H i column density A^h i = 
2 X 10 20 cm -2 , a temperature T = 10 4 K, a Doppler parameter 
b — 40 km s _1 , no dust, and a radial velocity of V exp = 300 
km s _1 . The Lya line equivalent width, given in units of the 
Doppler width Avq, increases from to oo from the top-left to 
the bottom-right panel. 



sorption is around x ~ x( V exp ), as this frequency corresponds to 
the line center frequency in the shell frame. The effect of dust is 
similar to Fig.[S]for spherical halos: it suppresses the emission 
peak and broadens the absorption, leading to an asymmetric 
absorption profile. The fraction of escaping photons far from 
the line center is equal to exp(-r fl ), as expected (see discussion 
section l3.2.2> . 

To illustrate the variety of emergent line profiles from 
a dust-free expanding shell when varying the intrinsic line 



strength, we show in Fig. [2U|a family of emergent spectra with 
intermediate cases between intrinsic pure continuum (cf. Fig. 
I19> and intrinsic monochromatic emission (cf. Fig. I12> . For a 
sufficiently large Lya equivalent width W(Lya), the second red 
peak (feature la), the small excess at x = (feature 3) and the 
blue bump (referred to as feature 2) become clearly visible. 

6. Lya line profile morphologies - models and 
qualitative comparison with observations 

For a better overview over the different simulations presented 
here and the resulting variety of Lya line profiles, and for a first 
qualitative comparison with observations we present a sum- 
mary in Tabled 

Schematically we may classify the emergent Lya profiles 
and considered geometries in the three following main groups, 
1) static media and symmetric profiles, 2) expanding/infalling 
halos, and 3) expanding shells, which we shall discuss now. 

6.1. Static media, symmetric profiles 

For simple static geometries with an embedded source emitting 
a symmetric spectrum around x = (i.e. monochromatic line 
radiation, a symmetric line centered on the Lya in the restframe 
of the medium, or a flat continuum) the emergent Lya line pro- 
file remains symmetric with two peaks, whose position is given 
by x p ~ ±(0.9 - 1) x (aro) 1 ^ 3 and little or no flux at line center 
(case 1 in TableHJ cf. Neufeld 1990, Hansen & Oh 2006). The 
detailed line shape (FWHM of each part, the possible presence 
of "bumps" etc.) depends on the properties of the scattering 
medium in ways quantified in detail by Neufeld (1990) and ex- 
amined also by Richling (2003) for clumpy media. Breaking 
the symmetry of the source spectrum or introducing velocity 
shifts between the source and the scattering medium lead to 
asymmetric profiles, as well known (cf. Neufeld 1990). 

Emission peaks superposed to a continuum - e.g. from a 
source with a pure continuum or continuum plus line emission 
- are easily "destroyed" in presence of dust (case 2). In this 
case the resulting profile is close to a Voigt profile. However, 
due to radiation transfer effects and the assumed geometry and 
"aperture", an integrated line profile fit in this manner will un- 
derestimate the true value of the column density. 

Purely or close to symmetric Lya line profiles are quite 
rarely observed. E.g. the spectra of the Lya blob LAB2 at 
z ~ 3.09 of Steidel et al. (2000) observed by Wilman et al. 
(2005), Lya emitters around a radio galaxy (Venemans et al. 
2005), or two LBG in the FORS Deep Field (Tapken 2005) 
come close to this, and may hence be related to static (or 
close to) media. For LAB2, however, other explanations have 
been put forward by Wilman et al. (2005) and Dijkstra et al. 
(2005b) including outflows + "absorbing screens" and inter- 
galactic medium (IGM) inflow. 

Some nearby starbursts such as I Zw 18 and SBS 0335-052 
show bf "Voigt-like" profiles (cf. Kunth et al. 1994, Thuan & 
Izotov 1997, Mas-Hesse et al. 2003). If related to a continuous 
source embedded in an H i cloud (case 2), such profiles can 
be reproduced with sufficient dust. In this case the H i column 
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Table 1. Schematic overview over the different simulations, the resulting line profile morphologies, and observed objects showing 
qualitatively similar profiles. Column 1 describes the adopted geometry, col. 2 describes the source position and spectrum, 
col. 3 the predicted line profile morphology, col. 4 the corresponding figure in the paper, and col. 5 possible objects showing 
qualitatively similar profiles. 



Case 


Geometry 


Source 


Line profile 


Figure 


Observations 


1 


static medium 


embedded source, 
monochromatic 


2 symmetrical peaks 
atx p = ±(ar ) 1/3 


Fig.[I] 


LAB2? 1 


2 


static medium 


embedded source, 
continuum 


2 symmetrical peaks if no dust 
DLA profile not related to r if dust 


Fig.uni 


IZwl8,SBS0335-052? 2 


3 


static medium 


external source, 
continuum 


DLA profile 

+ faint diffuse component 


Fig.CU 


DLA 3 


4 


expanding/infalling halo 


central source, 
monochromatic 


asymmetric emission peak 


Figs.|5l|9l 




5 


expanding/infalling halo 


central source, 
continuum 


P-Cygni without dust 

asymmetric absorption profile if dust 






6 


expanding shell 


central source, 
monochromatic 


1-2 red peaks (main peak at ~ 2V* max ), 
one blue peak 


FigsEmTunni 


LBGs 4 , low-z starbursts 2 


7 


expanding shell 


central source, 
continuum 


P-Cygni without dust 

asymmetric absorption profile if dust 


Figs-Eona 


IZwl8,SBS0335-052? 2 



References: 1 Wilman et al. (2005), 2 Mas-Hesse et al. (2003), 3 Adelberger et aL (2005), 4 see e.g. Shapley et al. (2003) and Noll et al. (2004) 



density derived from simple Voigt profile fitting would under- 
estimate the true value. Other configurations can, however, also 
lead to the same profiles. To clarify such ambiguities, a detailed 
analysis taking all known constraints into account will be nec- 
essary. 

A static medium illuminated by a distant background con- 
tinuum source (case 3) produces well known Voigt absorp- 
tion line profiles, e.g. seen as damped Lya absorbers (DLA). 
In such a geometry Lya radiation transfer effects need not be 
considered. A case of "Lya fluorescence" of QSO radiation 
on a nearby absorbing system has been observed recently by 
Adelberger et al. (2005): they detected a double-peaked emis- 
sion superimposed on a DLA profile. This can be related to our 
diffuse spectrum illustrated in Figll II 

6.2. Expanding/infalling halos 

As already discussed by Zheng & Miralda-Escude (2002) and 
Dijkstra et al. (2005a) for monochromatic sources emitting at 
line center, such geometries give rise to asymmetric profiles 
with a redshifted (blueshifted) main peak for expanding (in- 
falling) matter, and eventually secondary features (case 4). The 
position of the main peak shows the same dependence as for a 
Hubble outflow (Loeb & Rybicky 1999): it moves away from 
line center with increasing Ah and decreasing velocity V max - In 
the case of uniform emissivity, the positions of the secondary 
peak and the point of minimal flux indicate the external veloc- 
ity (see Fig. 8 in Dijkstra et al. 2005a). 

For sources emitting a pure continuum, the Lya radiation 
transfer leads to normal (inverted) P-Cygni profiles for expand- 
ing (infalling) halos. Again, as for static media, the presence of 
dust is able to destroy the emission peak leading then to absorp- 
tion line profiles with more or less pronounced asymmetries 
(case 5). 

Dijkstra et al. (2005b) have argued for possible infalling 
Lya halos in the case of the Lya blob LAB2 of Steidel et al. 



(2000), and of a z = 5.78 galaxy discovered by Bunker et 
al. (2003). Numerous LBGs and LAEs show asymmetric red- 
shifted Lya profiles, which could in principle be related to ex- 
panding halos. To distinguish this geometry from expanding 
shells or other geometries additional observational constraints 
are needed. For z ~ 3 LBGs, for example, expanding shells 
provide a good description, as we'll discuss now. 

6.3. Expanding shells 

Emergent Lya profiles from an expanding shell can be rather 
complex: from P-Cygni profiles or asymmetric absorption pro- 
files for an intrinsic continuum spectrum (case 7) to emission 
line profiles made of one (or two) red peaks and possibly a 
fainter blue bump for an intrinsically monochromatic (line) 
spectrum (case 6). In all cases the profiles show an asymmetric 
red wing related to Lyo- scattering in the outflowing medium. 
When Lya is in emission, the velocity shift of the main peak 
of the red part of the profile corresponds to velocities between 
V exp and ~ 2V exp , where V e xp is the expansion velocity of the 
shell. For column densities Ah ~ 10 20 cirr 2 the main red emis- 
sion peak measures quite accurately twice the shell velocity. 

Interestingly, spectroscopic observations of Lyman break 
galaxies (LBG) at z ~ 3 agree very well with the expectations 
for an expanding shell geometry. Indeed the composite spectra 
of the large LBG samples of Shapley et al. (2003) show not 
only clear signatures of outflows, as testified by a shift between 
Lya emission and interstellar absorption lines (with a mean 
value of Av em _ a b s ~ 650 km s ). The low-ionisation interstel- 
lar absorption lines are found blueshifted by Av = -150 ± 60 
km s _1 with respect to the stellar, systemic redshift, whereas 
the Lya emission is found redshifted by Av = +360 km s _1 
(Shapley et al. 2003), i.e. quite precisely at twice the velocity 
measured by the interstellar (IS) absorption lines. As the latter 
are formed by simple absorption processes, i.e. by interven- 
ing gas on the line of sight between the continuous source (the 
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starburst) and the observer, whereas Lya is prone to complex 
radiation transfer effects, the above empirical finding together 
with our modeling insight indicates very clearly that on aver- 
age the distribution and kinematics of the ISM in these LBGs is 
well described by an expanding shell. This is quite naturally the 
most simple geometry to simultaneously explain IS absorption 
lines blueshifted by an outflow velocity V exp and an asymmetric 
Lya line redshifted by twice this velocity ! 

As these two spectral features are formed on opposite hemi- 
spheres, IS absorption lines on the side approaching the ob- 
server and Lya from photons scattered back from the receding 
shell (cf. Pettini et al. 2001), this implies that overall (or "on 
average" over an entire LBG sample) the ISM shell structure 
must be fairly spherically symmetric rather than e.g. strongly 
bipolar. Such a structure may indeed be expected in case of 
strong starbursts triggering efficient large scale outflows (see 
e.g. simulations of Mori et al. 2002). 

For a shell geometry the finding of a velocity shift of 
~ 3 V^xp between Lya and the IS absorption lines also requires 
a large enough H i column density (cf. above). This is fully 
compatible with the typical/average value of A^h of the z ~ 3 
LBG sample of (Shapley et al. 2003) expected from the stan- 
dard correlation between Nh and the extinction Eb-v- 

Finally we note that the relatively smooth, singly peaked 
Lya line profiles observed by Shapley et al. (2003) and oth- 
ers are not in contradiction with our predictions. As already 
explained above, there are several ways to avoid the multiple 
peaks predicted for certain conditions on the red side of Lya 
(e.g. Figs. |Hto[l6}: among them are low expansion velocities, 
high column densities, or large Doppler parameters and com- 
binations thereof. Actually, for the relatively low outflow ve- 
locity of V exp ~ 150 + 60 km s _1 deduced from the Shapley et 
al. (2003) sample, and taking into account their relatively low 
spectral resolution, no distinguishable multiple red peaks are 
expected. At least for this sample, the problem posed by Ahn 
(2004) to avoid multiply peaked Lya profiles is ill posed. The 
possible but generally minor blue peak, which is predicted but 
not observed, does not pose a problem; in these high redshift 
objects it may be too faint to be detected or actually suppressed 
by the intervening interstellar H i clouds in the galaxy. 

Local starbursts (see e.g. Kunth et al. 1994, 1998, Gonzalez 
Delgado et al. 1998, Mass-Hesse et al. 2003) present P-Cygni 
profiles with a faint secondary peak on the blue side super- 
imposed to the absorption feature (e.g. IRAS 0833+6517), 
or more or less symmetric absorption profiles (e.g. IZw 18, 
SBS 0335-052). Mass-Hesse et al. (2003) interpret the vari- 
ety of observed Lya profiles as phases in the time evolution 
of an expanding supershell. Few galaxies of those showing P- 
Cygni profiles have measurements of the relative velocities be- 
tween the interstellar absorption lines, Lya, and the ionised gas 
(traced e.g. by Ha). Although, e.g. in Haro 2 the shift between 
Lya and the interstellar lines is ~ 2V exp , their spectra are a pri- 
ori not in contradiction with our simulations, as velocity shifts 
< 3 V exp are obtained for neutral column densities in the shell 
smaller than 10 20 crrT 2 , in agreement with the value derived 
by Lequeux et al. (1995). Furthermore, the secondary peak on 
the blue side of the P-cygni profiles is easily reproduced by 
Lya radiation transfer through a superbubble. This is maybe 



the strongest indication for such a geometry. Detailed model- 
ing of these nearby galaxies, taking all available observational 
constraints and spatial information into account, is necessary to 
shed more light onto these questions. 

6.4. General comments 

As the above shows, the expected Lya profiles show already 
quite a complexity depending on the values of the main param- 
eters considered (geometry, H i distribution and column den- 
sity, velocity field, Doppler parameter, dust to gas ratio, and 
input spectrum) and important degeneracies are found. Overall 
it is quite clear from this study and earlier ones that radiation 
transfer effects are important to shape the emergent Lya profile. 
The interpretation of observed P-Cygni profiles or asymmetric 
emission lines, for example, requires full radiation transfer cal- 
culations rather than superpositions of gaussian emission and 
Voigt-like absorption profiles as undertaken in various stud- 
ies (e.g. Mas-Hesse et al. 2003, Wilman et al. 2005). Such a 
detailed comparison with observations will be presented else- 
where. 

Obviously the geometries adopted in this paper are strongly 
idealised when compared to real galaxies. Asymmetries, non- 
unity filling or covering factors, arbitrary geometries, and ef- 
fects of different viewing and opening angles have to be con- 
sidered for more realistic simulations. Although our code can 
treat arbitrary density distributions, we have so far only inves- 
tigated homogeneous media. For complementary approaches 
treating clumpy media see e.g. Neufeld (1991), and Hansen 
& Oh (2006). Furthermore, in addition to the simple limiting 
cases of input spectra (monochromatic line plus continuum), 
other cases such as intrinsic Gaussian emission lines must be 
accounted for. All of these effects can easily be modeled with 
our generalised 3D radiation transfer code, and will be dis- 
cussed in future applications. 

Beyond the Lya radiation transfer modeled here for simpli- 
fied "galaxian" geometries it should be recalled that Lya line 
profiles predicted in this manner are still a priori prone to alter- 
ations due to the transmission through the intervening IGM., As 
well known (e.g. Madau 1995) this effect can significantly alter 
observed Lya line profiles — especially on the blue side — for 
objects of sufficiently high redshift (typically z <: 3). Effects 
of the IGM have e.g. been discussed by Haiman (2004), Santos 
(2004), and Dijkstra et al. (2005b), and are not treated here. 
For most purposes the IGM effects can be computed without 
the need for true Lya radiation transfer calculations. 

7. Summary and conclusion 

We have developed a new general 3D Lya radiation transfer 
code allowing for arbitrary hydrogen density, ionisation & tem- 
perature structures, dust distributions, and velocity fields using 
a Monte Carlo technique (Sect. 0. Currently the code works 
on a cartesian grid and it can handle Lya transfer problems 
in objects with arbitrary hydrogen column densities. The main 
radiation-matter interaction processes included in the code are 
Lya scattering, dust scattering, and dust absorption. UV/Lya 
photon sources with an arbitrary spatial distribution and arbi- 
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trary spectra can be treated. The main direct observables pre- 
dicted by our simulations are "global" (i.e. integrated over all 
viewing angles) Lya line profiles, spectra for different viewing 
directions and opening angles, and surface brightness maps at 
different wavelengths (or integrated over Lya). 

The code has been tested successfully and compared 
against analytical results and results from other simulations 
(e.g. Neufeld 1990, Ahn et al. 2001, Zheng & Miralda-Escude 
2002, Richling et al. 2003, Dijkstra & Haiman 2005a) for var- 
ious geometrical configurations including plane parallel slabs, 
disks, and expanding or infalling halos, and for computations 
with and without dust. Overall an excellent agreement is found, 
except in cases where assumptions made in the analytical cal- 
culation break down (see Sect. QJ- 

With the aim of understanding and interpreting ultimately 
the variety of Lya profiles observed in nearby and distant 
starbursts, we have examined the Lyo- line profiles predicted 
for several simple geometrical configurations and their depen- 
dence on the main parameters such as the H i column density, 
temperature, velocity, and dust content. We have considered 
slabs with internal sources, disks, expanding and infalling halos 
(Sect. [3J, externally illuminated slabs (Sect.|4}, and expanding 
shells (Sect. For the source spectrum we have considered 
the limiting cases of pure monochromatic line emission, a pure 
continuum, or intermediate cases with different line strengths. 
The two latter cases have rarely been discussed in the literature 
so far. 

Schematically, the following morphologies are found (cf. 
Table profiles reaching from doubly peaked symmetric 
emission to symmetric "Voigt-like" (absorption) profiles in 
static configurations with increasing dust content, and asym- 
metric red- (blue-) shifted emission lines with a blue (red) 
counterpart ranging from absorption to emission (with increas- 
ing line/continuum strength) in expanding (infalling) media. 
In principle, symmetric or nearly symmetric profiles are only 
obtained for static geometries or for small systematic veloci- 
ties. In practice, however, this may be altered by Lya scatter- 
ing/absorption in the intervening IGM, especially at high red- 
shift. 

Some specific results are worthwhile summarising here: 

- Pure Lya absorption line profiles may be observed in "inte- 
grated"/global spectra of objects with static geometries. In 
this case the naturally arising double Lya emission peaks 
have been suppressed by dust absorption. Standard Voigt 
profile fitting of such profiles will significantly underesti- 
mate the true hydrogen column density (see Sect.@J. 

- It should be noted that normal (inverted) P-Cygni like Lya 
profiles can be obtained in expanding (infalling) media 
from objects without any intrinsic Lya emission, as a nat- 
ural consequence of radiation transfer effects redistribut- 
ing UV (continuum) photons from the Lya line center to 
the red (blue) wing (e.g. Figs. l9l l20l >. Adding dust in such 
cases will progressively transform the P-Cygni profile into 
a broad asymmetric absorption profile. 

- Our 3D simulations of expanding shells lead us to revise 
the earlier results of Ahn and collaborators (Ahn et al. 
2003, Ahn 2004) and to clarify the radiation transfer pro- 



cesses explaining the detailed shape of Lya profiles pre- 
dicted from expanding shells. In particular we found that, 
for sufficiently large column densities (Nh ^ 10 20 cm" 2 ), 
the position of the main Lya emission peak is quite gen- 
erally redshifted by approximately twice the expansion ve- 
locity (Sect. |3- This is in excellent agreement with the ob- 
servations of z ~ 3 LBGs from Shapley et al. (2003), which 
show that Lya is redshifted by ~ 2V exp , where V exp is the 
expansion velocity measured from the interstellar absorp- 
tion lines blueshifted with respect to the stellar (i.e. galax- 
ian) redshift. This finding indicates also that large scale, 
fairly symmetric shell structures must be a good descrip- 
tion for the outflows in LBGs. 

As already clear from earlier investigations, radiation trans- 
fer effects even in the simple configurations just mentioned give 
rise to a complex morphology of Lya line profiles and some 
degeneraciesare found e.g. between the main parameters such 
as the hydrogen column density A^h, the Doppler parameter 
b, and also the intrinsic emission line width. However, codes 
as the one developed and presented here allow one to under- 
take detailed Lya line profile fitting to determine the physical 
properties of the emitting gas and of the surrounding H i and 
dust and to quantify possible degeneracies. With such tools it 
is now possible to fully exploit the information encoded in Lya 
profiles without recurring to physically unmotivated superpo- 
sitions of emission and Voigt absorption line profiles, which - 
except in specific geometric configurations - generally do not 
properly account for radiation transfer effects. Such detailed 
studies of nearby and distant starforming galaxies will be pre- 
sented in subsequent publications. Various other applications, 
combining also our MCLya code with the 3D radiation trans- 
fer and photoionisation code CRASH (Maselli et al. 2003) are 
foreseen. 
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